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ABSTRACT 

Context. The rich open cluster M67 is known to have a chemical composition close to solar, and an age around 4 Gyr. It thus offers 
the opportunity to check our understanding of the physics and the evolution of solar-type stars in a cluster environment. 
Aims. We present the first spectroscopic study at high resolution, R ~ 50,000, of the potentially best solar twin, M67- 1 194, identified 
among solar-like stars in M67. 

Oh' Methods. G dwarfs in M67 (d ~ 900 pc) are relatively faint (V ~ 15), which makes detailed spectroscopic studies time-consuming. 

Based on a pre-selection of solar-twin candidates performed at medium resolution by Pasquini et al. (2008), we explore the chemical- 
00 , abundance similarities and differences between M67-1 194 and the Sun, using VLT/FLAMES-UVES. 

Working with a solar twin in the framework of a differential analysis, we minimize systematic model errors in the abundance analysis 
compared to previous studies which utilized more evolved stars to determine the metallicity of M67. The differential approach yields 
precise and accurate chemical abundances for M67, which enhances the possibility to use this object in studies of the potential 
peculiarity, or normality, of the Sun. 
. Results. We find M67-1194 to have stellar parameters indistinguishable from the solar values, with the exception of the overall 

L-J metallicity which is slightly super-solar ([Fe/H] = 0.023 ± 0.015). An age determination based on evolutionary tracks yields 4.2 ± 

1.6 Gyr. Most surprisingly, we find the chemical abundance pattern to closely resemble the solar one, in contrast to most known solar 
, twins in the solar neighbourhood. 

O ji Conclusions. We confirm the solar-twin nature of M67-1 194, the first solar twin known to belong to a stellar association. This fact 

allows us to put some constraints on the physical reasons for the seemingly systematic departure of M67-1194 and the Sun from 
most known solar twins regarding chemical composition. We find that radiative dust cleansing by nearby luminous stars may be the 
explanation for the peculiar composition of both the Sun and M67-1 194, but alternative explanations are also possible. The chemical 
similarity between the Sun and M67-1 194 also suggests that the Sun once formed in a cluster like M67. 

, 
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I — , 1 . Introduction 1997, Melendez et al. 2006, Takeda et al. 2007, Melendez et al. 

I/*} , , , 2009, Ramirez et al. 2009). Although these stars in general have 

^ . For at least four decades, searches have been conducted for fundamental parame ters very close to solar, recent advances in 
. ■ stars with properties very similar to the Sun (see Hardorp 1978, hi h _ acc differential abundance analyses have proven al- 
g Cayrel de Strobel et al. 1981; see also the review by Cayrel de mQSt ^ of them tQ haye chemical compos i tions s i igM y } but sys . 

O ' I ^;w d f/. eCent * SU ™ y 6 - g - Me u lendez * tematically, deviating from that of the Sun. 
Ramirez 2007). It would be important to find a star with physical 

. . | characteristics indistinguishable from those of the Sun, a "per- 
fect good solar twin", as defined by Cayrel de Strobel (1996). 
The reasons for this importance are both physical and technical. 



A special opportunity in the search for solar twins is offered 
by the old and rich open cluster M67. It has a chemical com- 
position similar to the Sun with [Fe/H] in the range —0.04 to 
<>l i ne reasons ror tms importance are Dotn pnysicai ana tecnmcai. „ „_ . . . , . ,-5 in 

j<] ■ t., • „ , ■ • r i j. ■ i j i . i 0.03 on the customary logarithmic scale normalised to the Sun 

r.^ Physically, the statistics of solar twins would certainly contribute /IT , , „ „. . /_ ' , _ . „„„„ , 

5-h j ,. r i n (Hobbs & Thorburn 1991, Tautvaisiene et al. 2000, Yong et al. 

; to our understanding of the uniqueness or normality of the Sun Randich et pace et pa et * 

(cf. Gustafsson 1998). Technically, a solar twin would be useful 



Its age is also comparable to that of the Sun: 3.5-4.8 Gyr (e.g. 
Yadav et al. 2008). M67 is relatively nearby (~900pc) and is 
only little affected by interstellar extinction, which allows de- 
tailed spectroscopic studies of its main-sequence stars. The de- 
pleted Li abundance of the Sun seems rather representative of 
moved before an analysis of the spectroscopic features of the r , . . , _ . . . ,, m , „„,„, 

i. j ir i. c j solar-twins in the Galaxtic held (Baumann et al. 2010). M67 



in setting zero points in the calibration of effective-temperature 
scales, based on stellar colours. Another use would be in the cali- 
bration of night-time reflectance spectroscopy of solar-system 
bodies, where the spectral component of the Sun must be re- 



body itself can be performed 

Numerous searches and accurate analyses have resulted in a 
small sample of solar-twin candidates (Porto de Mello & da Silva 



also seems to contain Li-depleted G stars (Pasquini et al. 1997). 
M67 thus offers good possibilities of finding solar-twin candi- 
dates for further exploration. Pasquini et al. (2008) (followed by 
~W offprint requests to: Anna Onehag, e-mail: a paper of Biazzo et al. 2009) recently searched the cluster for 
Anna . OnehagS f ysast .uu.se solar analogs, and listed ten promising candidates. 

* Based on data obtained at ESO-VLT, Paranal Observatory, Chile, Here we present an analysis of M67-1194 (NGC 2682 YBP 

Program ID 082.D-0726(A) 1 194, ES 4063, ES IV-63, FBC 2867, MMJ 5357, SAND 770, 
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Fig. 1. Left panel: Cross-dispersion cut showing the eight-fibre set-up. Right panel: The pipeline-reduced data provided by ESO 
(red) compared to our work. Below the two spectra a difference spectrum is displayed for a direct comparison of the two reductions 
(data taken 2009-04-03: OB ID 331412). 



2MASS J085 10080+1 148527), a cluster solar-twin candidate 
suggested by Pasquini et al. (2008). Our analysis is based on 
high-resolution observations with relatively high signal-to-noise 
(S/N) ratio. In Section [2] we discuss the observations and some 
aspects of the data reduction. Section [3] describes the analysis 
method and the determination of fundamental parameters (T e g, 
logg, [Fe/H] and In Section |4] we present the results of a 
detailed analysis of a number of chemical elements, and these 
results are compared to those obtained for known twins in the 
Galactic field. In Section[5] we present a new age determination 
for M67, and in Section|6]we discuss the results. 

2. Observations 

The observations of M67-1 194 were carried out with the multi- 
object spectrometer FLAMES-UVES at ESO-VLT in Service 
Mode in the spring of 2009 during a period of three months (18th 
of January - 3rd of April). The observations analysed here are 
part of a larger project with a main goal to study atomic diffusion 
in stars of M67 (082.D-0726(A)). In each observing block of the 
project, one fibre of the spectrograph system was positioned on 
M67-1194 in order to collect as many observations as possible 
of this faint G dwarf. 

We obtained altogether 18 h in 13 observing nights (23 in- 
dividual observations). The spectrometer setting (RED580) was 
chosen to yield a resolution of R = A/AA =47,000 (1" fibre) 
and a wavelength coverage of 480-670 nm. The typical signal- 
to-noise ratio (S/N) per frame is 36 pixel" 1 (as measured in the 
line-free region between 6426 A and 6430 A). Barycentric radial 
velocities were determined from the individual spectra yielding a 
radial velocity of — 33.8±0.4km/s. This is in exellent agreement 
(la) with the mean radial velocity of the cluster as determined 
by Pasquini et al. (2008), v Iad = -33.30 ± 0.40, and Yadav et 
al. (2008), w rad = -33.673 ± 0.208. The latter study classifies 
M67-1194 as a high-probability proper-motion cluster member 
(99%). The frames were subsequently co-added for highest pos- 
sible S/N ratio. 

2. 1. Data reduction 

While ESO provides pipeline-reduced data for an initial assess- 
ment of the observed spectra, these are not intended for a full sci- 
entific analysis. In fact, we found that the extraction of our spec- 
tra by the ESO pipeline was problematic, mainly because the in- 
strumental setup of FLAMES-UVES combined with our choice 
of targets is challenging. Having eight fibres tightly packed in 
a slit-like fashion, FLAMES-UVES can obtain echelle spectra 
of up to 8 targets in parallel. However, the spacing between the 
individual fibres is rather small, which leads to non-negligible 



light contamination between adjacent fibres. This was made even 
more challenging by the fact that our target is considerably 
fainter than the one of the adjacent fibre. The ESO pipeline had 
severe problems extracting such a spectrum, which is apparent 
from a strong and periodic pattern of low-sensitivity spikes. This 
is illustrated in Fig.Q] where we show a cross-order profile and 
a section of the pipeline-reduced spectrum of M67-1 194 (in red) 
from frame OB 331412 (2009-04-03) — an illustrative frame in 
this context. 

Given the problems in the pipeline-reduced data, we adapted 
the code of the echelle-reduction package REDUCE (Piskunov 
& Valenti 2002) to perform an accurate spectrum extraction. 
Here we used special calibration frames provided by ESO, the 
so-called even and odd flat-field frames. In these frames, only 
every second fibre is illuminated, which eliminates crosstalk 
between fibres and allows us to calculate an accurate two- 
dimensional model of the individual aperture shapes. The ex- 
traction of the target spectra was then performed in two steps. 
We first used a very narrow aperture to extract the central (non- 
overlapping) part of each fibre. This allowed us to estimate the 
relative contribution of each individual fibre at each wavelength, 
which we then combined with the earlier determined aperture 
shapes to create a full two-dimensional model (which includes 
spectral features) for each fibre in the observed frame. We then 
eliminated the fibre-to-fibre contamination by subtracting the 
models of adjacent apertures prior to a full extraction of each 
individual spectrum. We show an example of a spectrum ex- 
tracted with this method in Fig. Q] (black). The periodic pattern 
of spikes has been removed successfully, and the agreement with 
the pipeline-reduced spectra outside the spikes is excellent. 

Fig.|2]shows a portion of the co-added spectrum covering the 
Mg lb lines. The S/N in the co-added spectrum reaches a mean 
of ^Opixer 1 . 

3. Analysis 

We performed a line-by-line differential analysis of the com- 
bined spectrum of M67-1 194 relative to the Sun. The latter was 
represented by a day-time spectrum of the sky taken with the 
580 nm setting of FLAMES-UVES in 2004 (Program ID 60. A- 
9022). The S/N of the solar spectrum is generally higher than 
that of M67-1194 (typically 230pixel _1 ). Together with the 
Kitt-Peak Solar Atlas (Kurucz et al. 1984), this made it possi- 
ble to identify blends of weak lines, in the wings of the lines 
measured as well as their surrounding regions used for setting 
the continuum. 

For analysis we used SIU (Reetz 1991), a tool to visual- 
ize and compare observed and theoretical spectra. SIU has a 
built-in line-synthesis module which uses one-dimensional hy- 
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Fig. 2. Observations of the Mg lb triplet region, for both M67-1 194 (black) and the FLAMES-UVES Sun (red). A difference spec- 
trum is plotted below. 
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Fig. 3. The Ha line observed in the solar spectrum (Kitt Peak 
National Observatory, KPNO) and in the spectrum of M67-1 194. 
Synthetic spectra with effective temperatures offset by ± 200 K 
are shown as solid lines. The most prominent telluric lines are 
marked out as vertical lines. 
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Fig. 4. The line-by-line difference in Fe I abundance plotted as a 
function of the excitation energy of each line. A slope insignifi- 
cantly different from zero is found, suggesting an effective tem- 
perature close to 5780 K. The standard deviation is indicated by 
the grey-shaded area. 



drostatic model atmospheres in LTE with an ODF representation 
of line opacity (MAFAGS, Fuhrmann et al. 1997, Grupp 2004). 
The line-formation code also assumes LTE. The highly differen- 
tial character of the analysis make our results very little depen- 
dent on the model atmospheres used; if, e.g., the more detailed 
MARCS models (Gustafsson et al. 2008) were used, the same 
results would be obtained. 

In a first step, we analysed the solar spectrum. We selected 
iron lines from the work of Korn et al. (2003), complemented 
by lines used by Melendez et al. (2009, priv. comm.) to cover 
a wide range in excitation energies. The KPNO solar atlas was 
consulted when setting the local continuum around lines of in- 
terest. With T cS = 5777 K, log g = 4.44, [Fe/H] = 0.00, a solar 



microturbulence of £ t = 0.95 km/s was found to minimize trends 
of abundance with line strength for both Fe I and Fe II. 

The derivation of a differential elemental abundance is a 
three-step process. First, the spectral region around every spec- 
tral line of interest is carefully normalized by comparing the 
relative fluxes in the solar and the stellar spectrum. Next, the 
solar spectrum is analysed varying the product of the elemen- 
tal abundance and the gf value until a satisfactory fit to the 
line profile is achieved. For strong lines damping parameters are 
based on quantum-mechanical calculations and adopted from 
VALD (Kupka et al. 1999). The Mg \b triplet line, used in the 
determination of the surface gravity, is fitted to the solar line 
varying the damping. The treatment of Ha follows Korn et al. 
(2003). The "external" profile, including effects of macroturbu- 
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Fig. 5. The line-by-line difference in Fe I and Fe II abundances 
plotted as a function of the equivalent width of the lines. Fe I 
lines return a slightly positive dependence (solid line), Fe II 
(dashed line) a slightly negative one suggesting an average £ t 
very close to the solar value of 0.95 km/s. 



lence, projected rotational velocity and instrumental profile, is 
assumed to be Gaussian with a width determined in this fitting 
process. Finally, the stellar spectrum is analysed and the abun- 
dance difference noted. Steps 2 and 3 are not carried out one after 
the other, to avoid biases that interactive differential line fitting 
may be subject to. Rather, a set of lines is analysed in the solar 
spectrum, then the same set of lines in the stellar one, without 
preconceptions as regards the abundance. 

We performed test calculations to quantify to which extent 
this interactive procedure gives results in agreement with those 
of a fully automatic procedure. An unweighted x 2 on the whole 
line profile was found to give systematically somewhat higher 
abundances in the analysis of the stellar spectrum. This is not 
surprising, as the experienced spectroscopist will tend to com- 
pensate for suspected blends while x 2 will blindly fit them as 
part of the line profile. However, the offset in abundance never 
exceeded 0.01 dex. Given the good agreement between these two 
fitting techniques, we are confident that the results of this study 
are independent of the line-fitting method used. 

3. 1. Fundamental stellar parameters and chemical 
composition 

The stellar parameters of M67-1 194 have been constrained em- 
ploying several spectroscopic techniques. The primary method 
to estimate the effective temperature is based on the self- 
broadened wings of Ha. As Fig.[3]shows, the agreement between 
the Ha profile observed for M67-1 194 and that of the Kitt Peak 
atlas is excellent. Within the noise level, there are no appreciable 
departures, except for the few telluric lines present in our obser- 
vations. The estimated effective temperature of M67-1 194 when 
based on the wings of Ha is 5780 ± 50 K. 

We furthermore analysed the excitation equilibrium of neu- 
tral iron (see Fig.[4]) and obtained a T e g equal to 5780, omitting 
one uncertain high-excitation line. We derived a mean metallic- 
ity for M67-1194 of [Fe/H]! = 0.023 ± 0.015 (la). The uncer- 
tainty in [Fe/H] translates into an uncertainty in effective tem- 
perature at the 20 K level (la). 
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Fig. 6. Observations of the gravity-sensitive wing of 
5183.6 A shown together with synthetic spectra for models 
log g values of 4.24, 4.44, and 4.64 (cgs units). 
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Combining the above two indicators, we adopt 5780 ± 27 K 
as our current best estimate of the effective temperature of M67- 
1194. The error is the combined error from the two methods 
propagated into the mean. 

Since the effective temperature is of key significance for the 
abundance results to be discussed below, we have performed two 
independent checks. One is to apply the line-ratio method of 
Gray (1994). Here, the central line depths of two neighbouring 
spectral lines, one at low and one at relatively high excitation, 
are compared. Following Gray, we have chosen the V I 625 1 .8 A 
line (x = 0.29 eV) and the Fel 6252.6 A line (x = 2.40 eV) 
which, assuming a solar [V/Fe] (consistent with the results of 
our analysis), suggests a temperature difference A = T c ff(M67- 
1194)-r Gff (0) = 60±70K. 

A second method is to use the observed colour of the star. 
According to Yadav et al. (2008), (B - V)ii 94 = 0.667 ± 0.01. 
Together with the comprehensive assessment of the reddening 
towards M67 (E(B - V) = 0.041 ± 0.004, Taylor 2007), this 
points towards a (B — V) of 0.626 ± 0.011 which is in the 
lower end of the interval (B — V)q — 0.62 — 0.65 suggested 
by contemporary studies (see Holmberg et al. 2006, who, how- 
ever, favour a high value around 0.64). If we adopt the recent 
calibration of (B — V)(T c g ) by Casagrande et al. (2010) we get 
T eff (M67-1194) = 5832 K. If we adopt a stellar K magnitude 
from the 2MASS catalogue we find (V — K) = 1.664, correct 
for reddening following Cox (2000) which suggests (V — K)q = 
1.556, and apply the Casagrande et al. (V — K) calibration we 
obtain T cff (M67-l 194) = 5773 K. Taking uncertainties in con- 
tinuum positioning for the line-depth estimates and in reddening 
for the colour estimates into consideration, all these results are 
consistent with our adopted value of 5780 K. The line-depth ra- 
tio method and photometry suggest that this value might rather 
be slightly underestimated than overestimated. 

The surface gravity is set by requiring lines of Fe II to return 
the same differential abundance as lines of Fe I (LTE ionization 
equilibrium). As in the case of Fe I, we opted for a careful selec- 
tion of relatively few well-observed weakly blended lines. Fig. [5] 
demonstrates the ionization equilibrium is established at log g 
= 4.44. Given a la scatter among the log g-sensitive Fe II lines 
(0.013 dex), the value of the surface gravity is uncertain at the 
0.04 dex level. In this step, the microturbulence is found to be 
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£t = 0.95 km/s (cf. Fig. |5), i.e. the same value as used in the 
analysis of the solar spectrum. 

Further evidence as regards the solar-like surface gravity is 
obtained from the analysis of the damping wings of the Mg lb 
lines. The Mg abundance was derived from the Mg I doublet 
lines at 6718.7 and 6319.2 A which are on the linear part of the 
curve of growth, as well as from the moderately strong Mg I 
5711 A line. These lines consistently give the result [Mg/Fe] 
= 0.02 ± 0.01. There are no other suitable atomic Mg lines in 
the spectral regions observed here. To get further constraints on 
the Mg abundance, we explored the wavelength region around 
5134-5137 A with Q 2 (23), Qi(23) and R 2 (ll) features of the 
(0,0) bands of the A-X system of MgH. This region is severely 
affected also by lines from the C2 Swan bands. The observed 
spectrum appears very similar but with about 8% stronger MgH 
bands, as compared with the solar spectrum. Our spectrum syn- 
thesis, assuming solar isotopic ratios, shows that the observed 
MgH bands are consistent with our abundance estimate for Mg 
and the somewhat enhanced overall metallicity. The adoption of 
a significantly lower Mg abundance (by 0. 1 dex) would, how- 
ever, lead to too weak MgH features. Adopting [Mg/H] = 0.02, 
the pressure-broadened red wing of Mg I 5183 A (the strongest 
line of the Mg I triplet with the largest log g sensitivity) indicates 
a log g of 4.44 ± 0.05 (see Fig. |6j. To this we add in quadrature 
the uncertainty stemming from the Mg abundance. This yields 
an error for log g of 0.06. We finally adopt 4.44 ± 0.035 as our 
current best estimate of the surface gravity of M67-1194. The 
error is the combined error from the two methods, Fe i/ll and 
Mg lb, have been propagated into the mean. 

Abundances are determined using the line list of Melendez 
et al. (2009), which was chosen in order to make our results di- 
rectly comparable to the twin results of that group. In the case 
of carbon, we needed to depart from that list replacing lines not 
available in our wavelength range with C2 lines. In the case of 
oxygen, inspection of the wavelength regions around the [Ol] 
lines at 6300.3 and 6363.7 A proved this region to be far too 
noisy for an accurate oxygen abundance determination. Within 
this noise, the spectrum does not depart significantly from the 
solar spectrum, and fits of synthetic spectra suggest the oxygen 
abundance ratio [O/Fe] to be solar with an estimated maximum 
error of about 0. 15 dex. The high excitation O I 6158. 1 line was 
found to be more useful at this level of S/N, allowing an esti- 
mated [O/Fe] = 0.07 ± 0.07. These errors are here to be seen 
as maximum errors, and the uncertainties are mainly related to 
the placement of the continuum. We note that this result is con- 
sistent with the mean [O/Fe] = 0.01 ± 0.03, as determined for 
ten somewhat brighter and hotter main-sequence stars in M67 
by Randich et al. (2006) using the 6300.3 A line. 

The resulting abundances are illustrated in Fig. [7] as differ- 
ences between the solar values and those of M67-1 194, and com- 
pared to the mean results for the solar twin sample of Melendez 
et al. (2009). We note that, given the similarity of the spectrum 
of M67-1194 to the solar one, uncertainties in line data effec- 
tively cancel. Similarly, errors due to inadequacies in the model 
atmospheres, e.g., in describing the formation of molecular lines 
in the upper layers of the atmospheres, should cancel, provided 
that the elemental abundances are practically the same, and no 
other effects but those related to the stellar fundamental param- 
eters characterize the atmospheres. 

4. M67-1194 in comparison with other twins 

The resulting abundance ratios for M67-1 194 relative to iron and 
normalized on the solar abundances, are plotted in Fig. UJ versus 
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Fig. 7. Abundance as a function of condensation temperature 
(according to Lodders 2003). Filled circles: this paper, plus 
signs: Melendez et al. 2009. Shaded area: this paper, dotted line: 
Melendez et al. 2009. The shaded area summarizes the result of 
this paper when submitted to a linear-regression analysis, giving 
A[X/Fe] = a+mT con d- The area delineates the range of regres- 
sion lines due to the uncertainties in a and m. The elements in 
the range 1300-1500K not identified are (from left to right): Cr, 
Si, Co, Ni and V. The error bars indicate the line-to-line scatter 
in abundances, e.g. the internal errors. If this error is below 0.02, 
we adopt 0.02 dex (indicated in the upper right corner and repre- 
sentative of all points without error bars). In the case of Na (two 
lines giving abundances in numerical agreement) we adopted an 
estimate of the fitting error. The effects of errors in fundamental 
stellar parameters lead to vertical shifts of the points by typically 
0.00-0.02 dex. 



the 50% condensation temperature of the various elements, ac- 
cording to Lodders (2003). Here, the results of Melendez et al. 
(2009) for their sample of eleven solar twins in the Galactic field 
are also indicated, and following those authors we plot the differ- 
ences between the solar abundances and corresponding results 
for the stars. 

Melendez et al. found that for almost all their stars, the 
volatile elements are relatively less abundant than in the Sun, 
while the refractory elements are somewhat enhanced in the 
stars. This effect is seemingly not present in the M67 twin, at 
least not to the same extent. It is thus seen (cf. Fig. [7]) that 
the abundance profile of M67-1194 is more similar to the Sun 
than for all the Melendez et al. twins, with one or two excep- 
tions. We perform a simple statistical test, essentially judging 
the probability that all black dots in Fig. [7] just by chance fall 
above the red dashed line for T con< \ > HOOK, and below it for 
T cond < 1400 K. We find that the probability for M67-1194 to 
be drawn from the same sample of stars as the majority of the 
Melendez et al. twins is less than 5%. On the other hand, the 
errors in the analysis are fully compatible with solar abundance 
ratios for the M67 star. This astonishing result will be discussed 
further below. 

One could ask to which extent the result shown in Fig. [7] 
would change if the fundamental parameters of M67-1194 
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Fig. 8. M67-1 194 compared to the FLAMES-UVES spectrum representing the Sun (red). The difference spectrum is shown below 
(offset by +0.7 units). At the current S/N level, there is no systematic difference between the two spectra at the position of the Li I 
doublet. 



are in error. We have found that the only error which could 
qualitatively affect the result would be an overestimate of the 
stellar temperature by 50 K or more. If the star is cooler by that 
amount, the abundances of the volatiles are overestimated so 
that when correcting for that error the points for C, O and S in 
the diagram would be pushed upwards towards the red dashed 
"Melendez twin" line while some points representing refracto- ; 
ries in the right part of the diagram would be pushed down. In ; 
view of the discussion of possible errors in T c fj(M67-l 194) we'- 
find, however, such an overestimate improbable. 
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4.1. Lithium content 

Pasquini et al. (2008) estimated the lithium content of ten 
solar-type stars in M67 including M67-1194 from FLAMES- 
MEDUSA data. The authors derived lithium abundances sim- 
ilar to that of the Sun from the Li doublet at 6707.8 A, using 
spectra with a resolution of R ~ 17,000 and a S/N ratio of 80- 
1 lOpixel" 1 . The weakness of the Li I 6707.8 A doublet in M67- 
1194 and the noise level of the spectra make it difficult to de- 
termine an accurate Li abundance for the twin. A major prob- 
lem is the location of the local continuum around the line. With 
our preferred positioning of the continuum (cf. Fig. [9] lower 
part) we find a stronger 6707 A line than the solar one, and a 
Li abundance of log e = 1 .26, corresponding to [Li/H] = +0.2 
dex. This continuum positioning presumes that the wide depres- 
sion around 6707.8 A is a real spectral feature, and not due to 
some instrumental effect. If the latter is assumed, a more local 
continuum drawing is appropriate (see upper part of Fig . [9j> . This 
would lead to a Li abundance close to solar, [Li/H] =0.0, i.e. a 
value closer but still greater than that obtained by Pasquini et al. 
(2008), [Li/H] < -0.2. We regard the Li I doublet line identifica- 
tion in M67-1 194 to be a probable but not a definitive one. It can 
therefore not be excluded that the true abundance is significantly 
below solar. Probably, however, the Li abundance of M67-1 194 
is solar or larger, but not exceeding the solar value by more than 
a factor of 2.5. 



T • 



0.98 - 



0.96 



z 0.94 



0.92 



0.90 




• M67-1 194 
o UVES-Sun 



6707.2 



6707.4 



6707.6 6707.1 
Wavelength [A] 



6708.0 



6708.2 



Fig. 9. Upper part: The spectrum of M67-1194 (black bullets) 
around Li I 6707 A with a local continuum compared to the 
FLAMES-UVES spectrum representing the Sun (red circles). 
The best theoretical fit to the spectrum is shown as a solid line. 
Lower part: A global continuum placement for M67-1194 to- 
gether with the best theoretical fit (both offset vertically by 0.045 
units). 



5. The age of M67 

Analysing a target very similar to the Sun, we profit from a num- 
ber of advantages. The spectroscopic stellar parameters are more 
accurately determined than is generally achieved. This is again 
a result of the strictly differential approach relative to the Sun, 
which is expected to promote cancellation of the modelling er- 
rors that usually dominate these types of analyses. The absolute 
accuracy in the values obtained is directly based on the accurate 
knowledge of the solar parameters. 
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This is then also the virtue of our age determinations. The 
stellar-evolutionary tracks are admittedly uncertain, due to e.g. 
diffusion, mixing and overshooting, the choice of mixing length 
parameter (a) and He mass fraction (Y), as well as the absolute 
values of the solar CNO abundances. However, a and Y are set 
to fit the solar radius and age. Assuming that a and Y are close 
to solar also for the solar twin, we may then use the tracks cali- 
brated in this way to obtain the age of the twin and thus of M67. 

Using the current version of the Victoria stellar evolution- 
ary code, we have computed several grids of tracks in order to 
evaluate the dependence of the predicted age on T e g, log g, Y 
and [Fe/H]. The main improvements that have been made to this 
code since the VandenBerg et al. (2007) study of M67 are the 
following ( 1 ) The gravitational settling of helium and turbulent 
mixing below envelope convection zones are now treated, using 
methods very similar to those employed by Proffitt & Michaud 
(199lfl (2) the latest rates for H-burning reactions (see Weiss 
2008, Marta et al. 2008) and the improved conductive opaci- 
ties given by Cassisi et al. (2007) have been adopted and (3) 
the assumed solar abundances are those reported by Asplund et 
al. (2009). To reproduce the properties of the Sun at its present 
age (4.57 Gyr; Bahcall et al. 2005), it is necessary to adopt 
an initial helium content corresponding to Y = 0.2553, along 
with a value of «mlt = 2.007 for the usual mixing-length 
parameter. The evolutionary tracks in each grid, for masses 
in the range 0.98 < M/M)q < 1.06, were evolved from 
fully-convective structures high on the Hayashi line to the so- 
lar age. Interpolations between the tracks yielded values of T c s 
and log g at ages ranging from 3.5 to 4.57 Gyr, which encom- 
passes most of the ages that have been determined for M 67 in 
the past decade (see below). An example of one of the computed 
grids on the (T c g, log <?)-plane is shown in Fig.fTOl 

From cubic spline interpolations in these results, the inferred 
age of the solar twin is 4.21 Gyr and its mass is estimated to be 
1.013 Mq. At T cff = 5780 K, one can readily determine that the 
gravity uncertainty implies a range in age from ss 2.6 to 5.6 Gyr. 
The age range spanned by the T c ff error bar is much smaller, 3.8- 
4.6 Gyr. Similar analyses of the model grids which were com- 
puted for slightly different helium and metal abundances yield 
the results that are listed in Table Q] The measured uncertainty 
in [Fe/H] (±0.015) implies an age uncertainty of approximately 
±0.2 Gyr, whereas the evolutionary calculations predict about a 
±0.3 Gyr change in age if the assumed value of Y is altered by 
±0.01. Melendez et al. (2009) propose that the Sun could have 
a reduced fraction of refractories in its present relatively thin 
convection zone, as compared with the composition of its inte- 
rior. The effects of such a chemical inhomogeneity not only for 
the Sun but also for M67-1194 on its age determination would 
cancel, but it is also possible that the twin, in spite of its so- 
lar abundance profile in the atmosphere, does not have a corre- 
sponding inhomogeneity. In that case, we find (using the results 
given in the bottom row of Table Q} a systematic error in the age 
determination, leading to an underestimate of the age by about 
0.35 Gyr. Table [2] provides a convenient summary of the errors 
in the fundamental parameters of M 67-1 194 and their effects on 
its inferred mass and age. 



At the solar age, settling reduces the surface helium abundance to 
Y — 0.2384, which is slightly less than the value inferred from helio- 
seismology, but similar to the values predicted by other Standard Solar 
Models, see the Bahcall et al. paper. Following Proffitt & Michaud, the 
free parameter in the adopted treatment of turbulent mixing was set so 
that our model for the Sun predicts the observed solar Li abundance. 



Our best estimate of the age of M 67-1 194 (4.2 ± 1.6 Gyr) 
agrees very well with with determinations based on alternative 
methods. Most fits of isochrones to the M67 color-magnitude 
diagram, whether using Garching (see Magic et al. 2010), 
Montreal (Michaud et al. 2004), Padova (Yadav et al. 2008), 
Teramo (Bellini et al. 2010), Victoria (VandenBerg et al. 2007), 
or Yale-Yonsei (Yadav et al.) isochrones, favour an age of 
4.0 ± 0.3 Gyr. Moreover, Bellini et al. have recently found that 
the age obtained from the cluster white dwarfs is fully consis- 
tent with the aforementioned turnoff age. However, higher ages 
cannot yet be completely excluded (see examples provided by 
Yadav et al. and Magic et al.). We note that the solar age is within 
the 2a error bar of M67's age. 

Table 1. The chemical-composition dependence of the inferred 
mass and age of M67-1 194. 



[Fe/H] 




Y 


M [Me] 


Age [Gyr] 


0.0 


0.01323 


0.2553 


1.002 


4.43 


0.023 


0.01413 


0.2450 


1.039 


3.87 


0.023 


0.01393 


0.2553 


1.014 


4.17 


0.023 


0.01375 


0.2650 


0.991 


4.44 


0.023 


0.01656" 


0.2553 


1.028 


3.50 


0.046 


0.01468 


0.2553 


1.027 


3.91 



Notes. A mixture with enhanced CNO abundances; specifically, 
[C/Fe] = [N/Fe] = [O/Fe] = 0.10. 



Table 2. Errors (e) in mass and age estimates for M67-1 194 due 
to uncertainties in fundamental stellar parameters. 



Fund, param. 


Value 


e{M/M Q ) 


e(Age) in Gyr 


T e s 


5780 ± 27 K 


±0.01 


T0.4 


logs' 


4.44 ±0.035 


±0.01 


Tl.5 


[Fc/H] 


0.023 ±0.015 


±0.01 


T0.2 


Y 


0.2553 ±0.013 


to.oi 


±0.3 



6. Discussion 

Our main finding, illustrated in Fig. [7] is that M67-1194 is in- 
deed more similar to the Sun than most of the nearby solar twins 
in the field. This is interesting in view of the possible expla- 
nations for the systematic departures of the field twin from the 
Sun in abundances with condensation temperature, found and 
discussed by Melendez et al. (2009, see also Gustafsson et al. 
2010). Basically, three possible explanations were discussed: (1) 
The solar proto-planetary gas disk was sufficiently long-lived, 
after the gas had been depleted in refractories by the formation 
of planetesimals and planets, for the solar convection zone to 
become so thin that the accreted gas could give visible traces 
in the solar spectrum. The field twins, however, in general ac- 
creted their gas disks much earlier. (2) The solar outer convec- 
tion zone never reached deep enough to encompass most of the 
solar mass, and was therefore much more easily polluted by the 
rarified infalling disk. This property of the Sun, presumably re- 
flecting the initial conditions of the cloud to later form the Sun 
and the Solar System or the details of the more or less episodic 
accretion history of the Sun (cf. Baraffe & Chabrier 2010), was 
not usual among solar-type stars and is thus not shared by most 
twins. (3) The proto-solar cloud was cleansed by radiation fields 
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Fig. 10. Plot of the 3.5 to 4.57 Gyr segments of the post-main- 
sequence tracks (dashed) that were computed for masses from 
0.98 to 1.06 A4q (indicated below each track) and the indicated 
values of Y and [Fe/H]. Open circles represent the T c gand log g 
values at the ages noted just to the right of the y-axis; and hence 
the solid curves are isochrones for those ages. The filled circle 
and attached error bars indicate the location of the solar twin, 
M 67-1 194. 



from massive stars, pushing dust grains out of the cloud. These 
effects must then have been greater for the solar cloud than for 
most of the twins. 

If the tendency found here for our M67 twin is representa- 
tive for the cluster, it speaks against the first explanation above. 
We see no reason why the gas disks around stars in a cluster 
should have longer life-times than those in the field. It might 
be that the initial conditions of star formation are different in a 
cluster than in the field, and that this might couple to the evo- 
lution of the stellar convection zones, but this is only specula- 
tion. Maybe, a more episodic accretion history, gradually build- 
ing stars from quite small entities, is more characteristic of clus- 
ter environment. In any case, there is independent evidence from 
radioisotopic abundances in the solar system that the Sun was 
born in a cluster (Looney et al. 2006, see also Portegies Zwart, 
2009) containing relatively close supernovae or massive AGB 
stars (Trigo-Rodriguez et al. 2009). 

The similarity of the age and overall composition of the Sun 
with the corresponding data of M67, and in particular the agree- 
ment of the detailed chemical composition of the Sun with that 
of M67-1194, could suggest that the Sun has formed in this 
very cluster. According to the numerical simulations by Hurley 
et al. (2005) the cluster has lost more than 80% of its stars by 
tidal interaction with the Galaxy, in particular when passing the 
Galactic plane, and the Sun might be one of those. We note that 
the orbit of the cluster encloses, within its apocentre and pericen- 
tre, the solar orbit. However, the cluster has an orbit extending 
to much higher Galactic latitudes, presently it is close to its ver- 
tical apex at z = 0.41 kpc (Davenport & Sandquist 2010), while 
the Sun does not reach beyond z = 80 pc (Innanen, Patrick & 
Duley 1978). Thus, in order for this hypothesis of an M67 ori- 
gin of the Sun to be valid, it must have been dispersed from the 
cluster into an orbit precisely in the plane of the Galactic disk, 
which seems improbable. 



7. Conclusions 

We draw the following conclusions from the analysis of M67- 
1 194 presented here: 

- Within the framework of a differential analysis, the stel- 
lar parameters (T c g, \ogg, mass and age) of M67-1194 are 
found to be compatible with the solar values, while [Fe/H] 
is found to be slightly super-solar. M67-1 194 is thus the first 
solar twin known to belong to a stellar association. 

- The chemical abundance pattern of M67-1194 closely re- 
sembles the solar one, in contrast to most known solar twins. 
This suggests similarities in the formation of M67-1 194 and 
the Sun. A common origin of both stars as members of M67 
is conceivable, albeit not likely, considering their different 
galactic orbits. If the chemical abundance pattern reflects en- 
vironmental effects, then the Sun was likely born in a cluster 
similar to M67. 

- Of the three scenarios for the peculiar abundance pattern of 
the Sun relative to field solar twins presented in Melendez 
et al. (2009), a particularly long-lived gas disk is the least 
probable due to the cluster membership of M67-1194. In 
the standard picture of stellar evolution with an early fully- 
convective phase, dust cleansing by nearby luminous stars 
seems to have affected both the Sun and M67-1 194. Another 
possibility is that the initial boundary conditions or accretion 
history of the early Sun were similar to those in the young 
cluster, and that these led to a shallow convection zone in 
the early evolution of these stars. More stars in M67 must be 
investigated to give further support to such explanations. 

Methodologically, we consider our method of age determination, 
based on spectroscopy of solar-like stars, to be an interesting al- 
ternative to other determinations of cluster ages. Its main advan- 
tages is that the uncertainties in cluster distance and reddening 
are not important. The age errors, in particular due to errors in 
the surface gravity, could be further reduced by higher S/N, a 
wider spectral range and observations of more stars. The uncer- 
tainty due to the unknown He abundance could be reduced by 
accurate studies of the horizontal branch or of visual binaries. 
In fact, the method can be extended to age estimates for other 
clusters, including globular ones. The basic idea is to make a 
systematic spectroscopic comparison with well known standard 
stars with very similar parameters in the nearby field or in other 
clusters. Even if the absolute ages of these standards are not as 
well established, the method offers an accurate age ranking, re- 
vealing possible age differences between different populations 
of clusters. 

In general, the highly differential method used here has great 
possibilities in other situations when detailed comparisons are to 
be carried out e.g. between Thick and Thin Galactic-disk stars, 
or different types of Halo stars. Then, a close adherence to a pair- 
wise selection of stars with similar effective temperatures and 
metallicities will make it possible to keep the systematic mod- 
elling errors at a low level. 

Acknowledgements. We thank the anonymous referee whose constructive criti- 
cism considerably improved the paper. AO acknowledges support by Varmlands 
nation (Uddeholms research scholarship). AJK acknowledges support by the 
Swedish Research Council and the Swedish National Space Board. 

References 

Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 
Bahcall, J. N., Basu, S., Pinsonneault, M., & Serenelli, A. M. 2005, ApJ, 618, 
1049 



8 



Onehag et al.: M67-1 194, an unusually Sun-like solar twin in M67 



Baraffe, I. & Chabrier, G. 2010, arXiv: 1008.4288vl 

Baumann, P., Ramirez, I., Melendez, J., Asplund, M., Lind, K. 2010, 
ai'Xiv: 1008.0575 

Biazzo K., Pasquini, L., Bonifacio, P., Randich, S. & Bedin, L. R. 2009, MmSAI, 
80, 125 

Bellini, A. et al. 2010, A&A, 513, A50 

Casagrande, L., Ramirez, I., Melendez, J., Bessell, M., Asplund, M. 2010. A&A, 
512, 54 

Cassisi, S., Potekhin, A., Pietrinferni, A., Catelan, M., & Salaris, M. 2007, ApJ, 
661, 1094 

Cayrel de Strobel, G., Knowles, N., Hernandez, G., & Bentolila, C. 1981, A&A, 
94, 1 

Cayrel de Strobel, G. 1996, A&ARv, 7, 243 

Cox, A. N. (ed.) 2000, Allen's Astrophysical Quantities, Fourth Ed., Springer- 

Verlag New York, Inc. 
Davenport, J. R. A. & Sandquist, E. L. 2010, ApJ, 71 1, 559 
Fuhrmann, K., Pfeiffer, M., Frank, C, Reetz, J., & Gehren, T. 1997, A&A, 323, 

909 

Gray, D. F. 1994, PASP, 106, 1248 
Grupp, F. 2004, A&A, 420, 289 
Gustafsson, B. 1998, SSrv, 85, 419 

Gustafsson, B., Edvardsson, B., Eriksson, K., J0rgensen, U. G., Nordlund, A., & 

Plez, B. 2008, A&A, 486, 951 
Gustafsson, B., Melendez, J., Asplund, M., Yong, D. 2010, Ap&SS, 328, 185 
Hardorp, J. 1978, A&A, 63, 383 
Hobbs, L. M. & Thorburn, J. A. 1991, AJ, 102, 1070 
Holmberg, J., Flynn, C, & Portinari, L. 2006, MNRAS, 367, 449 
Hurley J. R., Pols, O. R., Aarseth, S. J., & Tout, C. A. 2005, MNRAS, 363, 293 
Innanen, K. A., Patrick, A. T., Duley, W. W. 1978, Ap&SS, 57, 511 
Korn, A. J., Shi, J., & Gehren, T. 2003, A&A, 407, 691 

Kurucz, R. L., Furenlid, I., Brault, J., & Testerman, L. 1984, Solar Flux Atlas 

from 296 to 1300 nm, Kitt Peak National Solar Observatory 
Lodders, K. 2003, ApJ, 591, 1220 

Looney, L. W., Tobin, J. J., & Fields, B. D. 2006, 652, 1755 

Magic, Z., Serenelli, A. M., Weiss, A., & Chaboyer, B. 2010, ApJ, 718, 1378 

Marta, M., et al. 2008, Phys. Rev. C, 78, 022802 

Melendez, J., Dodds-Eden, K., & Robles, J. 2006, ApJ, 641, 133 

Melendez. J. & Ramirez, I. 2007, ApJ, 669, 89 

Melendez, J., Asplund, M., Gustafsson, B., & Yong, D. 2009, ApJ, 704, L66 
Michaud, G., Richard, O., Richer, J., & VandenBerg, D. A. 2004, ApJ, 606, 452 
Pace. G., Pasquini, L., & Francois, P. 2008, A&A, 489, 403 
Pasquini, L., Randich, S., & Pallavicini, R. 1997, A&A, 325, 535 
Pasquini, L. Biazzo, K., Bonifacio, P., Randich, S., & Bedin, L. R. 2008, A&A, 
489, 677 

Piskunov, N. E. & Valenti J. A. 2002, A&A 385, 1095 
Portegies Zwart, S. F. 2009, ApJ, 696, 13 

Porto de Mello, G. E, da Silva, R.. & da Silva, L. 2000, ASPC, 213, 73 
Porto de Mello, G. F. & da Silva, L. 1997, ApJ, 482, 89 
Ramirez, I., Melendez J., & Asplund, M. 2009, A&A, 508, 17 
Proffitt, C. R., & Michaud, G. 1991, ApJ, 371, 584 

Randich. S., Sestito, P., Primas, E, Pallavicini, R., & Pasquini, L. 2006, A&A, 
450, 557 

Reetz, J. K. 1991, Diploma Thesis, Universitat Miinchen 

Takeda, Y, Kawanomoto, S., Honda. S., Ando, H., & Sakurai, T. 2007, A&A, 

468, 663 
Taylor, B. J. 2007, AJ, 133, 370 

Tautvaisiene, G.. Edvardsson, B., Tuominen, I., & Ilyin, I. 2000, A&A, 360, 499 
Trigo-Rodrlguez, J. M., Garcla-Hernandez, D. A., Lugaro, M., Karakas, A. I., van 

Raai, M., Garcia Lario, P., & Manchado, A. 2009, M&PS, 44, 627 
VandenBerg, D. A., Gustafsson, B., Edvardsson B., Eriksson, K. & Ferguson, J. 

2007, ApJ, 666, 105 
Weiss, A. 2008, Physica Scripta, T133, 014025 

Yadav, R. K. S, Bedin, L. R., Piotto, G., Anderson, J., Cassisi, S., Villanova, S., 

Platais, I., Pasquini, L., Momany, Y. & Sagar, R. 2008, A&A, 484, 609 
Yong D., Carney, B. W., & Teixera de Almeida, M. L. 2005, AJ, 130, 597 



Appendix A: Chemical elements 



Table A.l. Fe I and Fe II lines used in the analysis. 



Fei 


[Fe/H] 


Fen 
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S414 073 


OA 

U.UH- 


S99S S^S 


u.uz 


<\4?^ 9S7 


09 
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U.UJ 


S9QS 9Q0 
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0^ 
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JJV I .tuu 


u.uu 






j j i y.j i h- 


u.uz. 






SS99 AAf\ 


n nzt 






5638.262 


-0.01 






5662.516 


0.03 






5679.023 


0.01 






5701.514 


0.02 






6065.482 


0.02 






6151.617 


0.01 






6200.313 


0.01 






6213.429 


0.04 






6229.225 


0.02 






6252.554 


0.04 






6421.349 


0.02 






6498.937 


0.01 
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Fig. A.l. Upper panel: Comparison of equivalent widths of M67- 
1194 and the FLAMES/UVES sky spectrum representing the 
Sun. ID polynomial fits for both Fel and Fell are shown as 
solid and dashed lines, respectively. Lower panel: The percent- 
age equivalent-width difference of Fe I and Fe II lines between 
M67-1194 and the FLAMES/UVES sky spectrum representing 
the Sun. Note that the abundances are not derived from equiva- 
lent widths, but from fits of synthetic spectra. 



Sun-like solar twin in M67 



Table A.2. Chemical elements (except Fe) and spectral lines in 
the analysis. 



Element 


A /£"1 m m s *1 tint, 

A(Sun - M67-1194) 


Element 


A(Sun - M67-1194) 


Ci&C 2 




VI 




5380.337 


0.03 


5670.853 


0.02 


5135.570 


0.01 


6039.722 


0.05 


5141.210 


0.00 


6081.441 


0.01 


5147.694 


0.015 


6090.214 


-0.01 


Oi 




6199.197 


-0.03 


6158.190 


-0.07 


Cri 




Nai 




5238.964 


0.00 


6154.225 


0.03 


5247.566 


-0.02 


6160.747 


0.03 


5272.007 


0.05 


Mgl 




5287.200 


0.04 


5711.088 


-0.01 


5296.691 


-0.03 


6318.710 6 


0.00 


5300.744 


0.00 


6319.240 


-0.04 


Mm 




All 




6013.513 


0.00 


6696.018 


0.00 


6016.673 


-0.01 


6698.667 


0.02 


6021.819 


0.01 


Si 1 




Co 1 




5488.983 


0.02 


5247.911 


0.02 


5665.554 


-0.02 


5301.039 


-0.015 


5690.425 


0.08 


5342.695 


0.02 


5701.104 


-0.02 


5530.774 


-0.01 


6125.021 


-0.02 


5647.230 


0.02 


6145.015 


0.05 


6189.000 


0.04 


6243.823 


-0.015 


Nil 




6244.476 


-0.03 


5589.358 


0.01 


6721.848 


0.01 


5643.078 


-0.01 


6741.630 


0.005 


6086.282 


0.02 


Si 




6108.116 


0.00 


6046.000 


0.005 


6130.135 


0.005 


6052.656 


0.02 


6204.604 


0.00 


6743.540 


0.01 


6223.984 


0.03 


6757.153 


-0.02 


6767.772 


-0.02 


Cai 




6772.315 


0.02 


5512.980 


-0.02 


Cui 




5590.114 


-0.02 


5218.197 


0.00 


5867.562 


0.05 


5220.066 


-0.005 


6166.439 


-0.02 


Yn 




6169.042 


-0.01 


4854.867 


-0.07 


Sen 




4883.685 


-0.01 


5657.896 


0.00 


4900.110 


-0.04 


5684.202 


-0.02 


5087.420 


0.04 


6245.637 


0.02 






Til 








5113.448 


0.02 






5219.700 


-0.01 






5490.150 


0.00 






5866.452 


0.01 






6126.217 


0.00 






6258.104 


-0.01 






6261.101 


0.05 







Notes. 1 ' For this line the KPNO atlas was used to represent the Sun as 
the FLAMES-UVES spectrum suffered from telluric absorption at this 
wavelength. 
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